Quantum Zeno Effect Explains Magnetic-Sensitive Radical-Ion-Pair Reactions 
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Chemical reactions involving radical-ion pairs are ubiquitous in biology, since not only are they 
at the basis of the photosynthetic reaction chain, but are also assumed to underlie the biochemical 
magnetic compass used by avian species for navigation. Recent experiments with magnetic-sensitive 
radical-ion pair reactions provided strong evidence for the radical-ion-pair magnetoreception mech- 
anism, verifying the expected magnetic sensitivities and chemical product yield changes. It is here 
shown that the theoretical description of radical-ion-pair reactions used since the 70's cannot explain 
the observed data, because it is based on phenomenological equations masking quantum coherence 
effects. The fundamental density matrix equation derived here from basic quantum measurement 
theory considerations naturally incorporates the quantum Zeno effect and readily explains recent 
experimental observations on low- and high-magnetic-field radical-ion-pair reactions. 

PACS numbers: 82.20.-w 



I. INTRODUCTION 

Recent experiments [l], Q on the excitonic energy 
transfer taking place in photosynthetic reactions have 
provided a tangible glimpse of quantum coherence effects 
being at play in biological systems. The possibility of 
quantum physics underpinning biological systems beyond 
the structural aspect has been entertained for a long time. 
It has also been clear 0, H| that in the decoherence-prone 
biological environment, some sort of protection of the 
quantum coherence must be at work before any effects 
related to the latter can surface. We will here show that 
a familiar biophysical system, namely radical-ion pairs 
and their reactions, exactly fulfills the aforementioned 
requirements and exhibits effects known from quantum 
physics experiments on well-isolated atomic systems and 
accounted for by quantum measurement theory. 

Radical-ion pairs, created by a charge transfer from 
a photo-excited donor-acceptor molecular dyad, are cen- 
tral in the reaction chain taking place in the photosyn- 
thetic reaction center [B|, 0] . The magnetic interactions 
0, [H of the two unpaired electrons in the radical-ion- 
pair with external magnetic fields and internal hyper- 
fine magnetic fields add another layer of complexity in 
the charge-transfer chain-reactions that convert the ab- 
sorbed photon energy to chemical energy vital for further 
biological function. Besides photosynthesis, magnetic- 
sensitive recombination reactions of radical-ion pairs are 
also assumed to underlie the biochemical magnetic com- 
pass 0, used by avian and possibly other species for 
navigation in the geomagnetic field [ll|, [H, OHII- Fur- 
thermore, radical-ion-pairs are also understood to partic- 
ipate in charge transfer reactions in DNA helices [TH, [l(| . 
It is thus clear that radical-ion pairs are found in several 
systems of biological significance, and the complexity of 
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their dynamics has been attacked from a wide range of 
theoretical and experimental disciplines. 

The time evolution of radical-ion-pair reactions is gov- 
erned by two distinct processes. The singlet (unpaired 
electron spins anti-aligned) and triplet (unpaired electron 
spins aligned) states of the radical- ion-pair and the inter- 
conversion between them brought about by the magnetic 
interactions Q are one piece of the dynamics. The other 
is the spin-state dependent charge-recombination of the 
radical-ion-pair. Singlet pairs recombine only to singlet 
neutral molecules and triplet pairs recombine only to 
triplet neutral molecules. 

However, since the 1970's, the theoretical treatment 
(reviewed in [13]) of radical-ion-pair reactions has been 
based on semi-classical/phenomenological equations in- 
volving the density matrix describing the pair's quantum 
state, not unlike the early rate-equation description of 
matter-light interactions. The latter is now known to 
mask important effects related to quantum coherence, 
which become transparent only in light of the full-blown 
quantum-mechanical description of atom-photon inter- 
actions [l||. It is the recombination process that has 
so far been treated phenomenologically in a single den- 
sity matrix equation simultaneously describing the mag- 
netic interactions within the radical-ion-pair. In light of 
quantum measurement theory applied to radical-ion-pair 
recombination reactions, it is here shown that this bi- 
ological system manifests quantum coherence, quantum 
jumps, the quantum Zeno effect and in principle quan- 
tum correlations, that is, the full machinery of quantum 
information science effects and concepts. 

The article is organized as follows. In the following sec- 
tion we will present the shortcomings of the phenomeno- 
logical master equation describing radical-ion-pair reac- 
tions so far while in Section III we put forward the 
quantum-mechanically consistent master equation. In 
Section IV we comment on spin relaxation, the pres- 
ence of which has in many cases conspired with the phe- 
nomenological treatment of the dynamics to produce a 
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FIG. 1: A simple schematic of radical-ion-pair reaction dy- 
namics. A photon excites a donor-acceptor molecular dyad, 
which after a charge transfer process creates a radical-ion- 
pair. The magnetic interactions of the two unpaired electrons 
(two dots) with the external magnetic field and the molecule's 
magnetic nuclei bring about a coherent singlet-triplet conver- 
sion. This would go on forever if it were not for the charge- 
recombination process that takes radical-ion pairs to neutral 
products. Angular momentum conservation enforces spin se- 
lectivity of the recombination, i.e. singlet (triplet) radical-ion 
pairs recombine to singlet (triplet) chemical products. 



description consistent with experiments. In Section V 
we use the fundamental quantum dynamic description 
of radical-ion-pair reactions put forward in Section III 
to reproduce recent experimental data by Hore and co- 
workers that clearly manifest quantum coherence 
effects not accounted for by the phenomenological the- 
ory. Finally, based on quantum measurement theory, the 
derivation of the fundamental master equation describ- 
ing radical-ion-pair reactions, which is done along the 
lines presented [2(1 HH, HH, [23[ by Milburn, Wiseman and 
co-workers on a rather similar system, namely coupled 
quantum dots, is presented in Section VI. 



II. PHENOMENOLOGICAL MASTER 
EQUATION 

The phenomenological or semi-classical density matrix 
equation that has been used so far to describe the dynam- 
ics of radical-ion-pair reactions, schematically depicted in 
Fig. CQ is 



dp 
~dt 



The first term describes the unitary evolution due to the 
magnetic Hamiltonian TL m ag, and the other two terms 
take into account the population loss due to singlet and 
triplet state recombination, respectively. The dimen- 
sion of the density matrix p is determined by the num- 
ber and the nuclear spin of the magnetic nuclei in the 
radical-ion-pair. It is given by dim = 4n, where 4 is 
the spin multiplicity of the two unpaired electrons and 
n = (2ii + 1) x (2I 2 + 1) x ... x (2Ij + 1) is the spin mul- 
tiplicity of the j magnetic nuclei existing in the donor 
and acceptor molecules. The Hermitian operators Qs 
and Qt are the singlet and triplet projection operators, 
respectively, given by 



Qs = 1/4 -si -s 2 
Qt = 3/4 + si • s 2 



(2) 
(3) 



where Si and s 2 are the spins of the two unpaired elec- 
trons (obviously Qs and Qt have to be properly writ- 
ten as 4n-dimensional matrices operating in the 4n- 
dimensional Hilbert space of the particular radical-ion- 
pair in consideration). The singlet and triplet projection 
operators are orthogonal, QsQt = QtQs = and add 
up to unity, Qs + Qt = 1, where 1 is the 4n-dimcnsional 
unit matrix. Finally, being projection operators, Qs and 
Qt satisfy the relations Q| = Qs and Q T — Qt, respec- 
tively. 

It is usually considered that the radical-ion-pair starts 
out in the singlet state, so that we can write for the ini- 
tial density matrix p(t = 0) = Qsl n i where n is the 
nuclear spin multiplicity. We let S(t) = Tr{p(t)Qs} and 
T(t) = Tr{p(t)QT} denote the time-dependent probabil- 
ity to find the radical-ion-pair in the singlet and triplet 
state. Due to the structure of |T]), the trace of the den- 
sity matrix decays exponentially to zero, since from (|TJ) 
it follows that 



dt 



-2k s S - 2k T T 



(4) 



This is the main problem of the phenomenological equa- 
tion m, i-c in order to describe population loss from the 
radical-ion-pair space into the neutral chemical products 
space due to charge recombination, the normalization of 
the density matrix is forced to an exponential decay. All 
coherences and populations are consequently also forced 
"by hand" to follow an exponential decay, thus eliminat- 
ing the actual presence of quantum coherence effects. 

What has not been realized so far, however, it that 
the tunneling taking place in the charge-recombination 
process is fundamentally a selective quantum measure- 
ment continuously interrogating the radical-ion-pair's 
spin state. If treated as such, the recombination pro- 
cess reveals well-known effects from the realm of quantum 
measurements, such as the quantum Zeno effect 0, 25 1, 
long-lived quantum coherences and quantum jumps 26 1. 



-i[Hmag, p]-k S (pQ S+Q SP)-kT(pQT+Qrp) (1) 



Quantum Zeno effects appear in several physical sys- 
tems, some of which are very similar to radical-ion-pairs, 
like the ortho-para conversion in molecular spin isomers 
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27 1 , ultra-cold atom tunneling through optical potentials 
28j , or the suppression of transverse spin- relaxation due 
to spin-exchange collisions in dense alkali-metal vapors 
[29l l3dj|. In the latter case, atomic spin-exchange colli- 
sions, of the form Si • S2, where Si and S2 are the electron 
spins of the two colliding atoms, probe the atomic spin 
state. When the collision rate (measurement rate) ex- 
ceeds the intrinsic frequency scale of the system, which 
is the Larmor frequency of spin precession in the applied 
magnetic field, the effective decay rate of the spin coher- 
ence is suppressed, a phenomenon that has led to the de- 
velopment of new ultra-sensitive atomic magnetometers 
[3lL \32\ . The radical-ion-pair tunneling into the neutral 
state is essentially a scattering process [33|, 0] , not un- 
like atomic collisions, that performs a measurement of the 
pair's spin state, since tunneling to singlet (triplet) prod- 
ucts can only proceed if the two unpaired electrons are 
in the singlet (triplet) spin state. Quantum Zeno effects 
have been extensively analyzed in the literature [35l . l36l ] , 
both with respect to pertaining physical systems, as well 
as the general conditions leading to the quantum Zeno 
effect or its inverse, the anti-Zeno effect 



III. QUANTUM MEASUREMENT THEORY 
DESCRIPTION OF RADICAL-ION-PAIR 
REACTIONS 

Using standard quantum measurement theory [3!| it 
readily follows (as shown in Section VI in detail) that the 
fundamental density matrix equation describing radical- 
ion-pair reactions and the associated magnetic-field ef- 
fects is 



dp 
dt 



■i[Hmag, p] ~ k S {pQs + QsP ~ ^Q S PQ s) 

-MpQt + QtP-ZQtpQt) (5) 



which is the same as the (p}, apart from the terms 
ZQspQs and 2QxpQr- It is these terms that are re- 
sponsible for the quantum effects that become important 
in the parameter regime where the recombination rates 
are larger than the magnetic interactions frequency scale. 
Since Qs + Qx = 1 , © can be simplified to 

^ = -i[H ma g, p] ~ HpQs + QsP - ZQspQs) (6) 

where k = kg + kx- The physical interpretation of ([5]) 
is that both singlet and triplet recombination channels 
essentially "measure" the same observable, Qs, with a 
total measurement rate ks + kx- In other words, the 
singlet (triplet) recombination channel is interrogating 
the radical-ion-pair whether it is in the singlet (triplet) 
state. But being in the triplet is not being in the singlet 
(Qs + Qt = 1), thus both channels measure the same 
observable. This physical interpretation of the actual 
dynamics is completely lacking from the phenomenolog- 
ical equation (fT]). The density matrix equation ^ has 



the property that the normalization of the density ma- 
trix does not change, i.e. Tr{p} = S(t) + T(t) = 1 at all 
times. It is here noted that as early as 1976, Haberkorn 
[38| arrived at ([5]) based on semi-quantitative arguments, 
but did not further consider it, exactly because it does 
not seem to describe population loss from the radical- 
ion-pair space into the reaction products space due to 
recombination, since Tr{p} = 1 at all times. All works 
henceforth have used the phenomenological equation |TJ . 

To take recombination into account, we need another 
stochastic equation involving quantum jumps 39] out of 
the radical-ion-pair space into the neutral chemical prod- 
ucts space. These are given by the probability of the 
singlet and triplet channel recombination, 



PS 
Pt 



2k s (Qs)dt 
2k T (QT)dt 



(7) 



The physical interpretation of the above equations is the 
following: ps (j>t) is the probability that the radical-ion- 
pair recombines into the singlet (triplet) channel in the 
time interval between t and t+dt. Both ps and px are ob- 
viously positive-definite numbers (the expectation values 
of Qs and Qx range between and 1). In other words, 
Haberkorn obtained the correct density matrix evolution 
equation, missing, however, the quantum jump equations 
([7]) and their physical interpretation, as the conceptual 
underpinnings of quantum measurement theory and open 
quantum systems were lacking at the time. In Section V 
on the explanation of recent experimental data we will 
explicitly explain how the quantum dynamic evolution 
described by ^ and (0 is actually simulated. 

Thus the fundamental quantum dynamical evolution of 
radical-ion-pair recombination reactions is given by ([6]) 
together with l[7|). To summarize the fundamental dif- 
ference of the quantum description based on ^ and (JT]) 
versus the phenomenological (fT]) we note that (fT]) is such 
that it continuously removes population from the singlet- 
triplet subspace into the chemical product space, and 
does so independently for each recombination channel. 
In reality however, the actual evolution of the radical- 
ion-pair's spin state is given by ([6]) (both channels affect 
the spin state evolution) and the molecules recombine in 
random times and in jump-like fashion according to ijTj). 
Thus even if the radical-ion-pair starts out in the singlet 
state, a high triplet recombination rate will significantly 
affect the spin state of the molecule, well before it has 
time to acquire a triplet character through the magnetic 
interactions. This physical statement is not embodied in 
the phenomenological equation {TJ). 



IV. SPIN RELAXATION 

There are two points that need to be commented on be- 
fore proceeding to present numerical simulations and the 
explanation of recent experimental data. The creation 
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FIG. 2: Magnetic Field Effect, Experimental Data and Theoretical Calculation, (a) Measured [Ty| magnetic field effect for 
two different magnetic fields B — 49 fiT (red curves) and B = 39 p,T (blue curves), (b) Reproduction of data from quantum 
measurement theory based on (JSJ and Q with a simple one- nuclear-spin (spin- 1/2) radical-ion-pair model with a diagonal 
hyperfine tensor given by a x — 8 ^s" 1 , a y = 2 /is' 1 and a z = 0. The recombination rates used are ks = 0.05 /is -1 and 
kx = 3.5 ^is" 1 . (c) Reproduction of data from the phenomenological theory fT]) with the same parameters as in (b). (d) 
Measured p^jj ] transient absorption traces at zero and high magnetic field and (g) corresponding magnetic field effect, (e) and 

(h) Reproduction of data from quantum measurement theory based on (|6]) and ([TJ with the same parameters as in (b). (f) and 

(i) Reproduction of data from the phenomenological theory ([TJ) with the same parameters as in (b). For the reproduction of 
the high field data we had to consider the creation dynamics of radical- ion pairs (mentioned in Section IV) with k cr = 4 ^s" 1 
and a high field relaxation rate ks,sr = 1-0 /is -1 . As explained in Section IV, due to the spin relaxation present at high fields, 
the qualitative disagreement between data and phenomenological theory is not that pronounced as at low fields. 



of radical-ion pairs does not happen instantaneously, but 
at a given rate k cr . This can be modeled by adding a 
source term to the density matrix equations, k cr e~ kcrV , 
the time- integral of which is equal to 1. This is not a fun- 
damental point and will henceforth be neglected, i.e. we 
will assume that the rate k cr is much larger than all other 
rates of the problem, unless noted otherwise. Second and 
more important is the effect of spin-relaxation. This can 
be modeled by a term —k sr (p — po) , where k sr is the 
relaxation rate and pn = l/4n is the fully mixed density 
matrix, i.e. a diagonal 4n-dimensional matrix with di- 
agonal elements equal to l/4n. In the phenomenological 
equation ([T]) the analogous spin-relaxation term should 
read ~k sr p, since by design equation (fTJ) has p = as the 
infinite-time solution, hence the relaxation term should 
work towards the same steady state. We will "turn on" 
spin relaxation by setting k sr ^ 0, to show that in the 
presence of spin relaxation the results following from the 



quantum-mechanical description based on © and (JT]) 
naturally blend in with the ones derived from (TT]) . This is 
the main reason why so far the phenomenological density 
matrix equation {1} has provided a more or less consistent 
theoretical description of experimental observations. Ob- 
viously, the charge recombination process itself induces 
spin relaxation, and this is what is described by ([6]). By 
the rate k sr we describe all other spin relaxation effects, 
for example spin-lattice relaxation. It is also noted that 
diffusion or other effects add another layer of complexity 
to the dynamics described by ^ and ([7]) and will not 
be considered here. In other words, equations and 
([7]) are the fundamental dynamical equations describing 
magnetic-sensitive radical-ion-pair reactions in the ideal- 
ized case of negligible diffusion. This is not unrealistic, 
since at low temperatures (as for example in the exper- 
iment reported in Ref. [l!|) the effects of diffusion and 
more importantly, spin relaxation, are suppressed. 
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V. EXPLANATION OF EXPERIMENTAL DATA 

We will now show that the quantum-mechanical anal- 
ysis of the radical-ion-pair reactions put forward in Sec- 
tion III seamlessly explains recent experimental data, 
while the phenomenological theory is unable to do so. 
In the experiment recently reported [l9j by Hore and 
co-workers, a carotenoid-porphyrin-fullerene (CPF) triad 
was used to study the magnetic field effect on the radical- 
pair's recombination dynamics. We have performed nu- 
merical simulations using the simplest physically realiz- 
able radical-ion-pair model, that containing just one spin- 
1/2 magnetic nucleus. In this case the density matrix is 
8-dimensional and the magnetic Hamiltonian reads 

Umag = Uj(s lz + S 2z ) + I • A • Si (8) 

where we considered the external magnetic field to be 
along the z-axis, B = Bz, and the electron Larmor 
frequency parameter (the real Larmor frequency is two 
times higher) is to = j e B, with 7 = 1.4 MHz/G. The hy- 
perfine coupling tensor A couples the single nuclear spin 
existing in either the donor or the acceptor molecule with 
the corresponding unpaired electron. Anisotropic hyper- 
fine interactions embodied in the hyperfine tensor are 
responsible for the directional sensitivity of the reaction, 
i.e. the sensitivity of the reaction yields on the angle be- 
tween the external magnetic field and the local coordinate 
frame defined by the hyperfine tensor. In the experiment 
reported in fl~oj ] the angular dependence of the reaction 
has also been measured, and can also be reproduced the- 
oretically from the quantum dynamic description of the 
reaction, however, we will present this topic in a different 
manuscript. 

By direct numerical integration of either ([T|) or (JBJ) and 
([7]) we calculate the population of the triad molecules in 
the singlet-triplet subspace at time t. This population 
is what is measured in the transient absorption mea- 
surement. The " magnetic- field effect" is the difference 
(S(t) + T{t)) B - (S(t) + T(t)) B =o , i.e. the difference 
of the measured transient absorption at a magnetic field 
B from that at zero field. For the sake of completeness 
we note that while the phenomenological predictions are 
obtained by a direct integration of (|TJ) , the quantum dy- 
namic predictions based on ^ and the quantum jump 
equation ([7|) are obtained as follows: we start with Nq 
radical-ion pairs, and in each time step the number of 
existing radical- ion-pairs N(t) evolves according to 

N(t) = N(t-dt)(l-ps-p T ) (9) 

i.e. the number of radical-ion-pairs disappearing to the 
chemical product space (either singlet or triplet) between 
time t — dt and time t is N(t — dt) times the probability to 
jump to a recombined state (either singlet or triplet) dur- 
ing the time interval dt. This probability is just ps +Pt, 
where ps and px are given by ([7|) . Along the reaction, the 
density matrix p describing the spin state of the radical- 
ion-pair is evolved according to 



We take the radical-ion-pair to be in the singlet state 
initially, i.e. pit — 0) = Qs/2 (in general for a sin- 
glet initial state it is pit = 0) = Qs/n, where n is the 
nuclear spin multiplicity, which in this case of one spin- 
1/2 nucleus is n = 2). In Fig. 2a we show the measured 
magnetic field effect at low magnetic fields (earth's field) , 
while in Fig. [2)3 we fully reproduce the data using the 
quantum dynamics description of © and Q, whereas 
Figj2j: shows that the phenomenological equation {T]) is 
not capable of reproducing the observed magnetic field 
effect with the same parameters, and the same is true for 
any choice of the system's parameters. 

It is seen that the quantum dynamic description based 
on (|6|) and ((7|) agrees with the data both qualitatively 
and quantitatively (the magnitude of the effect), whereas 
the opposite is true for the phenomenological predictions 
shown in Fig. [2J;. It is noted that the time-scale of the 
measured magnetic field effect is missed by about a factor 
of 2, since in reality the molecule's dynamics are deter- 
mined by tens of nuclear spins, whereas here we have con- 
sidered just one nuclear spin. It is also important to note 
the following: in Fig. [2t> we reproduce the data using 
a small value for the singlet-channel recombination rate 
and a much larger value for the triplet-channel recombi- 
nation rate, while the molecule starts out in the singlet 
state. This is exactly the parameter regime where the 
quantum Zeno effect is actually manifested, i.e. as men- 
tioned before, the high triplet recombination rate has a 
tangible effect on the spin dynamics, as it appears in the 
total measurement rate k = ks+kx- This is a regime that 
the phenomenological theory fails to account for even at 
the qualitative level. 



We will now demonstrate that if we take into account 
a possible presence of spin-relaxation, the magnetic field 
effect is suppressed both in the phenomenological descrip- 
tion and in the quantum-mechanical, and both descrip- 
tions qualitatively merge into each other. The presence of 
spin-relaxation effects is one reason why the phenomeno- 
logical reaction dynamics of equation ([1]) have been rel- 
atively successful until now. This is expected since in 
general relaxation effects tend to equalize level popula- 
tions and damp quantum coherences, adversely affecting 
the precision of any spectroscopic measurement, such as 
the determination of the applied magnetic field. We will 
further show that the bi-phasic response observed in the 
data (Fig. [2^) and reproduced theoretically (Fig. [2)3) 
fades away when we "turn on" spin relaxation, and re- 
sembles the response produced by the phenomenological 
equation ([1]). The reason that the experiment by Maeda 
et al [lil manifests effects not accountable by the density 
matrix equation used until now is the suppressed pres- 
ence of spin-relaxation effects at the low temperatures at 
which the experiment was carried out. Indeed, in Fig. 
[3] we turn on spin-relaxation and it is evident that the 
magnetic field effect derived from (J6j) and ([7|) resembles 
that derived from ([T]) as the spin relaxation rate k sr is 
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FIG. 3: Magnetic Field Effect and Spin Relaxation. The magnetic field effect is simulated at a magnetic field B = 49 /iT. 
Traces a,b,c have been produced with the quantum dynamics given by © and ([7)l. while traces d,e,f have been produced by 
the phenomenological equation (JTJ). The same calculations have been done for three different values of the relaxation rate 
k sr — 0, 1, 10 /is -1 . Plot (a) is the same as plot (b) of Fig. [2] i.e. the one reproducing the low-field data of [jg| . It is seen 
that by increasing the spin relaxation rate, the two theories produce qualitatively and quantitatively similar results, i.e. for 
increasing k sr the bi-phasic response of the magnetic-field-effect is suppressed and the magnitude of the effect is roughly the 
same in both theories. It is thus clear that the presence of spin relaxation has in many cases masked the insufficiency of the 
phenomenological description and made it appear consistent with experiments. The other way around, it is clear that the low 
temperatures used in the recent experiment by Hore and co-workers pl| has suppressed spin-relaxation effects to the extent 
that the true quantum nature of radical-ion-pair reactions could surface. 
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increased. 

Finally, also the high-field data of can be explained 
with the same magnetic Hamiltonian used to reproduce 
the low-field data shown in Fig. [2K To this end we 
take into account the creation of radical-ion pairs, which 
does not take place instantaneously, but at a given rate 
k cr , modelled by the term dN/dt = Nok cr e~ kcrt , where 
No is the total number of molecules. In Fig. [2Ji we 
show the experimental results for the transient absorp- 
tion at a magnetic field of B=0 and B=8 mT, with the 
quantum measurement calculation shown in Fig. [2^ and 
the phenomenological prediction in Fig. [5f. Similarly, 
Figures Ek,h,i show the data and theoretical predictions 
for the magnetic field effect, i.e. the difference of the 
traces shown in Figures [2ji,e and f. The initial negative 
response is due to non-zero spin-relaxation at high mag- 
netic fields, which we model as described previously and 
which has been mentioned to exist in fl9j | . It is interesting 
to note the fact that the data reproduction by quantum 
measurement theory is capable of picking up this small 
spin relaxation (k sr = 1 /^s" 1 ) existing at high magnetic 
fields. Taking into account the simplicity of the magnetic 
Hamiltonian used for the calculations, the agreement be- 
tween experimental data and the predictions of quantum 
measurement theory is excellent, as both the absolute 
value of the transient absorption as well as the strength 
of the magnetic field effect and the time scale of the phe- 
nomenon are perfectly reproduced. It is also interesting 
to note that at high magnetic fields the disagreement be- 
tween data and the predictions of the phenomenological 
theory is not that pronounced as at low fields. This is 
exactly due to the small spin relaxation present at high 
fields. 




FIG. 4: Generic level structure and recombination dynam- 
ics in a radical-ion-pair, taking place in four steps: 1, pho- 
toexcitation, 2, radical-ion-pair creation, 3, tunneling-induced 
quantum measurement of the pair's spin state and 4, final de- 
cay to the ground state. Only the singlet reservoir of vibra- 
tionally excited states (of the neutral molecule DA) is shown 
for simplicity. Another such reservoir exists for the triplet 
chemical products. 



VI. DERIVATION OF MASTER EQUATION 

In Figure 4 we depict in more detail the generic model 
for the radical-ion-pair creation and recombination dy- 
namics shown in Fig. 1. A donor-acceptor molecule DA 
is photo-excited (D*A) and a subsequent charge-transfer 
creates the radical-ion-pair (D + A~). The singlet and 
triplet states of the radical-ion-pair (^D+A", T D+A~) 
are split by internal magnetic interactions of the radical- 
ion-pair's two unpaired electrons with external magnetic 
fields and internal hyperfine couplings. The radical-ion- 
pair is initially created in the singlet state, which is not an 
eigenstate of the magnetic Hamiltonian, and therefore a 
singlet-triplet (S-T) coherent mixing commences. A tun- 
neling event into a vibrationally-excited state of the neu- 
tral recombined molecule DA, which quickly decays into 
the ground state, completes the photon-initiated radical- 
ion-pair reaction. 

As is well known [4(1 El[ , electron transfer in radical- 
ion-pair recombination reactions is fundamentally a 
quantum-mechanical tunneling process. We will show 
that this process constitutes a continuous quantum mea- 
surement of the pair's spin state. Like every quantum 



measurement, this one is no exception to the rule that 
measurements performed on a quantum system lead to 
decoherence [42[ . However, under appropriate conditions 
involving the measurement rate and the intrinsic (mag- 
netic) frequency scale of the radical-ion-pair, the quan- 
tum Zeno effect 0, IH[ appears and leads to physical 
consequences completely against the intuition that has 
developed over the years based on the phenomenological 
description of (JXJ) . 

In the following, we are going to capitalize on the re- 
markably strong analogy between radical-ion-pairs and 
yet another physical system, namely two coupled quan- 
tum dots [lo, [H, 0, Sl| being continuously interro- 
gated by a point contact. This analogy is schematically 
depicted in Fig. O We are going to identify the anal- 
ogous physical observables of the two systems and then 
derive the corresponding evolution equation describing 
radical-ion-pair reactions. The electron hopping between 
the two dots is the analog of the S-T coherent mixing 
taking place in the radical-ion-pair, whereas the mea- 
surement performed by the point contact corresponds to 
the spin-state-dependent electron tunneling into an ad- 
jacent excited state of the neutral molecule. Radical-ion 
pairs in the singlet (triplet) state will tunnel to a nearby 
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magnetic 
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FIG. 5: Analogy between coupled quantum dots and radical- 
ion-pair reactions, (a) Schematic of coupled quantum dots 
interrogated by a double point contact (PC), (b) Schematic 
of radical-ion-pair with the singlet and triplet recombination 
reservoirs. 



empty reservoir of singlet (triplet) vibrationally excited 
states, the annihilation and creation operators of which 
are denoted by a (b) and (r (w), respectively. The cor- 
responding Hamiltonian is 



(10) 



As explained before, the radical-ion-pair consists of two 
unpaired electrons, one at the site of the donor (D) and 
one at the site of the acceptor (A). We denote by c and 
c* the fermion annihilation and creation operators of the 
electron at the acceptor site. Obviously, if (c^c) = 1 we 
have a radical-ion-pair, if however, (c'c) = 0, we have 
a neutral recombined molecule. The Hamiltonian of the 
radical-ion-pair thus reads 



Hrip C^c(e ~t~ H-mag 



(11) 



where TL m ag is the magnetic Hamiltonian operating in 
the spin space of the radical-ion-pair (for a single nuclear 
spin it is given in |8]) and it is readily generalized for 
any number of nuclear spins by adding more hyperfine 
interaction terms) and e the energy of the radical-ion- 
pair with respect to the ground state neutral molecule. 
Lastly, the tunneling Hamiltonian reads 



(12) 



The interpretation of the first term is that if the radical- 
ion-pair is in the spin singlet state it can charge re- 
combine, reducing the radical-ion-pair occupation num- 
ber to zero and creating an occupied vibrationally ex- 
cited state of the singlet reservoir. Similarly for the 
triplet state in the second term. In summary, the total 
Hilbert space consists of the singlet and triplet reservoirs, 
the occupation of the radical-ion-pair's acceptor electron 
and the spin space of the radical-ion-pair. For example, 
in the case of singlet-channel recombination, the initial 
and final quantum states coupled by Ht would look like 



I**) = l{0},l)lx) and |*/) = |{n},0}| Xs ), where the first 
set of quantum numbers {n} is the excitation of the vi- 
brational reservoir, the second is the electron occupation 
number of the acceptor site, and \x) is the spin state of 
the radical-ion-pair (two electrons plus nuclei), while |x s ) 
is the electron spin singlet state of the neutral donor elec- 
trons with the nuclear spin state being the same in \x) 
and \xs)- The tunneling amplitudes are [40| Ti = Vifi 
(i = a,b), where Vi is the electronic matrix element and 
fi the vibrational overlap between the nuclear wavefunc- 
tions of the reservoir states \i) and the radical- ion-pair. 

A very similar set of Hamiltonians has already been 
treated by Milburn, Wiseman and co-workers 0, |2l|, [22j 
for the case of quantum dots, and we will here proceed 
along the same lines to derive the master equation for 
radical-ion pairs. We write the total Hamiltonian as H = 
Ho + V, where 



Hq = C f ce + Hres 

and the perturbation V is 

V = C^cHmaq + H T 



(13) 



(14) 



and move to the interaction picture defined by Ho, in 
which the perturbation is Vi = e lW °* Ve~ lW °* and it fol- 
lows that 



Vi(t) = c)cH 



■mag 



+ c.c. 



y^T b cb^QTe 



i(ui b — e)t 



(15) 



The proof of (|15p is given in the Appendix. We will now 
apply time-dependent perturbation theory up to 2 nd or- 
der. The density matrix of the total system (spin state 
of radical-ion-pair + electron + reservoir states) is de- 
noted by W. The density matrix of interest describ- 
ing the radical-ion-pair spin state will be p, and will be 
calculated by tracing out the other degrees of freedom 
from W. During the time of coherent singlet-triplet mix- 
ing, i.e. as long as the radical-ion-pair has not recom- 
bined, the reservoir states are empty, which means that 
the density matrix of the singlet and triplet reservoirs is 
Pres = J2a l°)aa(0| ® ^6 1 ) & fc, <0 1 . For the same reason, 
the density matrix for the occupation state of the accep- 
tor electron is p c = |1} CC (1| and the total density matrix 
is then given by W(t) — p{t) ® p res ® p c . Obviously 
p{t) = Tr a , b , c {W(t)}. 

Since the tunneling part of the perturbation Vi is linear 
in the reservoir operators a, a\b and b\ it easily follows 
that in the l st -ordcr of the perturbation expansion (h = 

1) 

W(t + dt) = W(t) 

-idt[Vi{t),W{t)\ 

-dt [ Gfti[V/(«),[Vj(ti),W(t)]] (16) 
Jo 
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the terms in Vi containing the tunneling amplitudes T a 
and Tf, do not contribute, because e.g. Tr Q f,{p res a} = 0. 
Since Tr c {ctcW roa . g p G } = (c t c}H mag = W mag) it follows 
that the l st -order term in p^|) results in 

(%) =-i[H mag ,p] (17) 

\ at j 1 st_ ordcr 

We will now proceed to calculate the integral appearing 
in the 2 nd -order term of (fTB|) . In this term we will neglect 
T~tmag with respect to Ht, under the assumption that 
the measurement dynamics embodied in the tunneling 
Hamiltonian (|12p dominates the magnetic interactions, 
as is actually the case in the regime where the quantum 
Zeno effect is manifested. This approximation is vali- 
dated by the fact that the resulting master equation has 
the ability to reproduce experimental data. We will show 
the calculation just for the a-reservoir states (the singlet 
reservoir) in (|15p , as the term involving the triplet reser- 
voir (6-states) is dealt with in exactly the same way. In 
the integral appearing in (jTHJ> , containing Vj to 2 nd -order, 
terms involving T a from Vi(t) and T a from Vj(ii) contain 
the operator c , which is zero. The only terms that sur- 
vive are the ones containing T a from the V/(t)-term and 
T* from the complex conjugate part of the V7(<i)-term. 
These terms lead (after taking into account the fact that 
cross-terms involving different reservoir states a and a' 
as well as cross terms involving both reservoirs do not 
contribute) to an expression of the form 

-dtJ2 \T a \ 2 [ dhe^-^-^ica^Qs, [JaQ s , W}] 
+ c.c. (18) 

Taking into account the fact that for the empty reservoir 
states it holds (a^a) = and (aa^) = 1, while for the 
acceptor electron it is (c'c) = 1 and (cc 1 ) = 0, it follows 
that when tracing out the a and c degrees of freedom 
from the double commutator in (|18p , we are left with the 
operator pQs — QspQs- The complex conjugate (c.c.) 
term will then offer Qsp—QspQs, thus in total we obtain 
the Lindblad form pQs + QsP ~ %QspQs- We finally 
consider the resonance condition uj a ss e and make the 
usual approximation J d^ie l ^ a_£ ^ t_tl ' w 2nS(uj a — e), 
leading to the desired result, which after doing the same 
calculation for the triplet reservoir states b, reads 

J = -i[H mag , P ] - k s V[Q s ]p - k T V[Q T ]p (19) 

where the super-operator T>[B] acts on the density matrix 
p according to 

V[B]p = B^Bp + pB^B -2BpB' ! (20) 

We have thus derived ([S]), with the recombination rates 
given by (we have re-introduced h for the consistency of 



units) 

k s = (2TT/H 2 )Y / \Ta\ 2 S(e~LJ a ) (21) 

a 

k T = (2Tr/h 2 )Y,\Tb\ 2 S{e~u; b ) (22) 

b 

The resulting master equation ([6]) is indeed of the Lind- 
blad form that follows from standard quantum measure- 
ment theory [3^, |42|, |45[ , when the measured observable 
is Qs and the measurement rate is k = ks + kx, lead- 
ing to the physical interpretation that the spin-dependent 
recombination process is a continuous quantum measure- 
ment of the radical ion pair's spin state. 



VII. CONCLUSIONS 

In summary, we have shown that there exists a bio- 
logical system, radical-ion pairs and their reactions, in 
which the full machinery of quantum measurement the- 
ory can be fruitfully applied, leading to an understanding 
of quantum effects visible in the laboratory but masked 
by the phenomenological description of radical-ion-pair 
reactions employed until now. A fundamental quantum 
phenomenon, the quantum Zeno effect, familiar from 
carefully prepared atomic physics experiments, is seen 
to participate in the complex realm of biologically sig- 
nificant chemical reactions, further opening up the way 
to the emerging field of quantum biology, i.e. the ex- 
ploration of the manifestations of quantum phenomena 
and the application of quantum information concepts in 
biological systems. 
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APPENDIX 
1. Proof of equation (|15l) 

We will show that e lectct ce -* ectct = ce~ tet . Since c is a 
fermion annihilation operator, it is c 2 = 0, thus e lec ct c = 
c. Furthermore, since (c ! c) n = c^c for any integer power 
n, it follows that e~ %tc ct — e _let c^c + cc\ leading to 
the desired result. In a similar fashion it can be shown 
that e in -r CS t e -iH rcs t _ a t e *"a* an( j similarly for the b- 

degrees of freedom, thus follows equation (fT~5|) . 
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